Vortex Pinball Under Crossed AC Drives in Superconductors with Periodic Pinning 

Arrays 
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Vortices driven with both a transverse and a longitudinal AC drive which are out of phase are 
shown to exhibit a novel commensuration-incommensuration effect when interacting with periodic 
substrates. For different AC driving parameters, the motion of the vortices forms commensurate 
orbits with the periodicity of the pinning array. When the commensurate orbits are present, there 
is a finite DC critical depinning threshold, while for the incommensurate phases the vortices are 
delocalized and the DC depinning threshold is absent. 
PACS numbers: 74.60.Ge,74.60.Jg 



A wide variety of dynamical systems can be modeled as 
a classical particle moving in a periodic potential. Exam- 
ples of such systems are Josephson-junction (JJ) arrays 
fj],|J, sliding charge density waves (CDW) ||||], atomic 
friction ||, vortices in superconductors moving over pe- 
riodic substrates |6|-|l2|, electrons at low magnetic fields 
in antidot arrays |J.y, as well as ions In. 51 and colloids 
ful in optical trap arrays. When an AC drive is super- 
imposed over the DC drive, resonances or phase-locking 
effects can occur which have been extensively studied in 
J J arrays || and CDWs j| . Signatures of the phase lock- 
ing include oscillations in the depinning threshold with 
increasing AC drive amplitude, and locking of the parti- 
cle velocity to a fixed value over a range of the DC drive. 
In superconductors with periodic substrates, phase lock- 
ing of driven vortices has been experimentally observed in 
samples containing ID periodic modulations Q . More re- 
cently, experiments, simulations, and theory have shown 
that in superconductors with periodic hole arrays, phase 
locking or Shapiro steps occur when combined DC and 
AC drives are applied. 

In all these phase locking systems the AC and DC 
drives are in the same direction and the motion can be 
considered effectively one-dimensional (ID). Different ef- 
fects can arise in a two-dimensional (2D) system. Re- 
cently it was shown that new kinds of phase-locking ef- 
fects occur when the AC drive is applied perpendicular to 
the DC drive for systems with periodic substrates jlj] . 

Here we study vortex motion in superconductors with 
periodic arrays of pinning sites and three combined driv- 
ing forces. A DC drive fdc is applied along the longitu- 
dinal direction, and two AC drives which are 90 degrees 
out of phase are added in the longitudinal and transverse 
directions. In the absence of pinning and at fdc — 0, the 
vortices move in a circular orbit with radius and eccen- 
tricity determined by the amplitudes and frequencies of 
the AC drives. We focus on a regime just above the first 
matching field B^, where each pin is occupied by one vor- 
tex and a small number of additional vortices move in the 
periodic potential created by the pinned vortices. The 
existence of such interstitial vortices above Ba, has been 



inferred from transport measurements and direct imag- 
ing in samples with periodic holes B, and in this same 
regime Shapiro steps have previously been observed ||. 

The system we propose has similarities to the "electron 
pinball" model studied by Weiss et al. jl3| and others jl4) , 
in which classical cyclotron electron motion is induced by 
a magnetic field in samples with periodic arrays of anti- 
dots. In these systems, peaks and dips appear in the 
magnetoresistance as a function of field. The features 
are believed to arise when classically moving electrons 
follow pinned or commensurate circular orbits enclosing 
integer numbers of dots, in which the electrons can travel 
without scattering off the dot potentials. At the incom- 
mensurate orbits the electrons are scattered and diffuse 
throughout the sample. There are important differences 
between the vortex " pinball" system and the electron sys- 
tems. The vortices interact with a long-range smooth 
or egg-carton square potential, rather than simple point 
scatters, and thus the moving vortices must follow square 
orbits (in the case of square pinning arrays) , rather than 
circular orbits. A particularly attractive feature of the 
vortex system is that the shape of the orbit can be care- 
fully controlled experimentally by changing the ratio of 
the AC amplitudes, phases or frequencies. This allows 
new kinds of anisotropic orbits and commensuration ef- 
fects to be produced, which are not obtainable in the 
electron pinball models. 

We consider a thin superconductor containing 2D vor- 
tices, which is the appropriate model for the recent ex- 
periments in superconductors with hole and dot arrays. 
The vortex- vortex interaction has the form of a logarith- 
mic potential, U v = —A v hi(r), with the energy normal- 
ization A v = $q/87tA. Here, $o is the flux quantum and 
A is the effective 2D penetration depth for a thin film 
superconductor. The overdamped normalized equation 
of motion for a single vortex i is 



dt 



(1) 



where the damping term r/ is the Bardeen-Stephen fric- 
tion. The force from vortex i on the other vortices is 
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(e) (f) 
FIG. 1. Vortex positions (dots) and trajectories (lines) for 
fdc = 0.0 and different isotropic AC amplitudes, (a) A = 0.16, 
(b) A = 0.21, (c) A = 0.3, (d) A = 0.355, (e) A = 0.39, and 
(f) A = 0.47. 

f™ = — Yjj^i ^iUv ( r ij)' The long range interaction is 
treated with a fast converging summation method [fl7|| . 
The pinning force f vp arises from the pinning sites which 
are modeled as short range attractive parabolic wells. 
The pinning sites are placed in a square array of side 
L, and each pin has a radius of r p = 0.15L, which is 
within the typical experimental ratios of r p /L = 0.14 to 
0.3. The DC driving term f c i c is applied along the sym- 
metry axis of the pinning array, in the x-direction. The 
AC driving term is ?ac — -Asin^^^x + B cos(uiBt)y, 
where we fix wa/wb = 1-0. For most of the results pre- 
sented here we consider the case of a system containing 
64 pins and a vortex filling fraction of B/B^ = 1 + 1/64; 
however, we have found the same results for higher filling 
fractions (such as B/B$ = 2, 7/4, 3/2, and 5/4) at which 
the interstitial vortices form a symmetric pattern so that 
interstitial vortex interactions cancel. The initial vortex 
position is found by annealing from a high temperature 



with no driving and cooling to T = 0. The DC drive is 
increased in increments of 0.0001 and the sample is held 
at each drive for 3 x 10 5 time steps to ensure a steady 
state; the DC depinning threshold is determined from the 
time averaged vortex velocities < V x >. 

We first consider the case of isotropic AC driv- 
ing with A/B = 1.0 and show that commensurate- 
incommensurate depinning transitions occur as a func- 
tion of increasing A even at fdc = 0. Example vortex tra- 
jectories are illustrated in Fig. 1. For A = 0.16 [Fig. 1(a)] 
the vortex is confined to move in a small circular orbit in 
the middle of the plaquette. For A = 0.21 [Fig. 1(b)] the 
AC amplitude is large enough that the vortex encircles 
n = 1 pin in an orbit that is slightly tilted due to the 
counter-clockwise vortex movement through the square 
potential. In Fig. 1(c), for A = 0.3, the orbit encircles 
n = 4 pins. In Fig. 1(d), for A = 0.355, the orbit ra- 
dius is somewhere between n — 4 and n = 9, and the 
vortex interacts strongly with the occupied pins, scat- 
tering off them. The vortex becomes delocalized, and 
diffuses throughout the sample, avoiding areas near pins 
due to the repulsion from the pinned vortices. The orbit 
switches intermittently between n = 4 and n = 9. In 
Fig. l(e,f) we show localized vortex orbits for A = 0.39 
and 0.47, with n = 9 and n = 16, respectively. We 
also find a stable orbit with n — 25 for A = 0.57. For 
0.42 < A < 0.45, as well as for 0.52 < A < 0.55, the 
vortex is delocalized and moves in a manner similar to 
that shown in Fig. 1(d). 

In general we find localized orbits at values of A for 
which the vortex trajectory encloses n = m 2 pins, where 
to is an integer. By comparison, in electrons in antidot 
lattices, Weiss et al. |l5| ] found commensurate or pinned 
orbits when the number of dots encircled was n = 1, 2, 
4, 9, 16, and 21. We do not observe any stable orbits for 
n = 2 and 21. In the electron systems jl^,|lj] these orbits 
correspond to states where a portion of the electron orbit 
closely approaches the dots, which are point scatters. In 
the vortex system, due to the long range vortex-vortex 
interactions, the mobile vortices cannot follow orbits that 
approach arbitrarily closely to the pinning sites. 

In Fig. 2 we show the DC depinning threshold Fd p /F^ p 
vs A for the isotropic case A = B. Here, F® p is the 
DC depinning force for A = 0. The commensurate or- 
bits are pinned at low fd c , producing a series of peaks 
in Fd p at the commensurate AC values. At the incom- 
mensurate phases the depinning threshold vanishes. For 
0.0 < A < 0.2, the vortices move in the n = intersti- 
tial orbit illustrated in Fig. 1(a). Fdc decreases as the 
AC amplitude increases until reaching a minimum value 
at A = 0.19, followed by a peak in F dc for A = 0.225, 
which corresponds to an n = 1 orbit [Fig. 1(b)]. The 
second peak at A = 0.31 corresponds to the n = 4 orbit 
[Fig. 1(c)]. For the delocalized incommensurate orbits 
[Fig. 1(d)] of 0.34 < A < 0.36, F dp vanishes. A finite F dp 
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FIG. 2. The DC depinning force F dp /F% p vs AC amplitude 
A for the isotropic case A = B. F$ p is the depinning force 
for zero AC drive. Inset: average DC velocity < V m > vs 
DC driving force /d c for an incommensurate orbit A — 0.43 
(dotted line) and a commensurate orbit A = 0.37 (solid line). 

is regained when n ~ 9 orbits appear for 0.36 < A < 0.42. 
As A continues to increase, Fdp is nonzero for A values at 
which the stable n = 16 and n = 25 orbits are observed, 
with Fdp = in portions of the regions between these val- 
ues of A. We have checked the oscillatory behavior for 
different values of L and find the same general behavior. 
Since the radius of the vortex orbit R ~ A/ to a, a similar 
series of peaks in Fdc occurs if A is fixed and 1/oja is 
varied. In the inset of Fig. 2 we show typical curves of 
the average DC velocity < V x > vs fdc for an unpinnned 
incommensurate orbit at A = 0.43 and a pinned com- 
mensurate orbit at A = 0.37, illustrating the well defined 
sharp depinning threshold for the pinned orbits. 

The maximum Fdp values in Fig. 2 do not decrease with 
increasing A, indicating that all of the commensurate 
orbits are pinned equally well. In addition the bound- 
aries between the pinned and unpinned regions are sharp. 
Thus the behavior of Fd p vs A clearly differs from that as- 
sociated with Shapiro steps [p]PJlO[| , where the depinning 
threshold oscillates with A according to a Bessel function 
Jo(A), the peaks in the depinning force are smooth, and 
the peak Fdp value gradually decreases with A. 

We next consider anisotropic or elliptical orbits for A ^ 
B. In Fig. 3 we show F dp /F* p vs B/A for fixed A = 0.225 
and varying B, where F^ p is the depinning force for B = 
0. Again we find strongly pinned orbits, indicated by 
peaks in Fdp] however, in this case the pinned orbits occur 
when the vortex orbit encircles n = m pins. If a higher 
value of A is chosen such that the orbit encircles two 
pinning sites in the transverse direction in one period, 
commensurate orbits that encircle n = 2m pins appear. 
The peaks in the depinning curve shown here are much 
more symmetric than those in Fig. 2. This may be due 




FIG. 3. The DC depinning force Fd p /F dp vs AC amplitude 
B/A for the anisotropic case where A = 0.225 and B is varied. 
F*l p is the depinning force for B — 0.0. 

to the fact that the elliptical orbits are less perturbed by 
the square pinning potential than the isotropic circular 
orbits. 

In Fig. 4 we present several anisotropic commensurate 
orbits at the depinning peaks shown in Fig. 3. Here the 
orbits for Fig. 4(a,b,d) encircle n = 1, 2, and 3 pins, 
respectively. This trend continues for the higher peaks. 
In Fig. 4(c) we also show a sliding orbit at B — 0.3 with 
fdc above depinning. 

We now discuss experimental systems in which these 
phases can be observed. For superconductors with peri- 
odic antidot arrays, the pin geometry should be chosen 
such that the pins have vortex saturation numbers of 
one, so that above B^ the additional vortices will sit in 
the interstitial regions. For pins with higher saturation 
numbers, commensuration effects should still be observ- 
able at higher matching fields when interstitial vortices 
start to appear. Samples with very low intrinsic pinning 
should be used to enhance the effect. Samples in which 
Shapiro steps for moving interstitial vortices have been 
observed would be ideal. As in the case of Shapiro steps 
JTo| , the variation of the depinning force versus AC drive 
amplitude should be most pronounced at filling fractions 
where the interstitial vortices form a symmetrical pattern 
and the interstitial vortex interactions effectively cancel, 
such as at B/B^ = 2.0, 1.75, 1.5, and 1.25. Using samples 
such those in Ref. J9| where Shapiro steps are observed, 
for a pinning lattice spacing of a — 2fim, uj — 2ttv, and 
v = AQMHz, commensuration effects should be observ- 
able with applied crossed AC currents from to 10 I c 
where I c is the critical current. Superconductors with 
rectangular and triangular (rather than square) pinning 
arrays should also exhibit these phenomena; however, 
the stable orbits would encircle different numbers of pins 
than the orbits described here. In superconductors with 
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(a) 



(b) 



(c) (d) 
FIG. 4. Vortex positions (dots) and trajectories (lines) for 
fixed A = 0.225 and (a) B = 0.225, (b) B = 0.285 and (d) 
B = 0.35. In (c) a sliding orbit is shown for B = 0.3 and 

fdc > Fdp. 

magnetic dot arrays, the dots can act as long range re- 
pulsive sites, so that pinned orbits of the type described 
here could be observed at vortex filling fractions such 
as B/Bcf, = 1 and 1/2. In addition, superconductors 
with 2D smoothly modulated surfaces should also pro- 
duce commensurate and incommensurate orbits. Other 
promising systems in which these phases could be ob- 
served include 2D JJ arrays, 2D atomic friction models, 
and colloids on optical pin-scapes [jl6|. A possible ap- 
plication of our results would be particle segregation in 
multi-species systems, such as colloids on periodic sub- 
strates, where the different species have different mobili- 
ties. Here it should be possible to tune the AC drive such 
that one species would be pinned and another depinned. 
The species could then be segregated with a DC drive. 

To summarize, we have numerically studied the motion 
of vortices interacting with a 2D periodic potential cre- 
ated by immobile vortices located at pinning sites placed 
in a square array. We apply two AC drives, perpendicu- 
lar to one another and out of phase by 90° , such that the 
vortices move in a circle in the absence of a substrate po- 
tential. For AC drives of equal amplitude, as a function 
of increasing AC amplitude we find a series of pinned or- 
bits enclosing n = m 2 pinning sites. Each of these orbits 
has a finite depinning threshold to an additional applied 
DC force. At AC amplitude values between these pinned 
orbits, the vortices are delocalized and diffuse through 
the sample, and the DC depinning threshold is zero. Ex- 
perimentally these states can be observed as a series of 



pinned and non-pinned regions as a function of AC am- 
plitude or frequency. For anisotropic AC drives, we find 
a series of asymmetric pinned orbits which enclose n = in 
pinning sites. We call our model the vortex pinball model 
in analogy to the electron pinball system for electron cy- 
clotron motion in anti-dot arrays. We also suggest other 
systems in which these phases can be observed experi- 
mentally. 
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